Introduction to the model-free control of microgrids 



Loi'c MICHEL 

GREI - Departement de genie electrique et genie informatique 
Universite du Quebec a Trois- Rivieres 
CP. 500, Trois-Rivieres, G9A 5H7, Canada, QC 



Abstract 

This letteiQ presents the application of the model- free control approach to 
the microgrid control. We show in simulation that the method allows to 
control, with a simple controller, voltage, current and power of inverter- 
based microgrids. 
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1 Introduction 



The model-free control methodology, originally proposed by [T] , has been widely suc- 
cessfully applied to many mechanical and electrical processes. The model- free control 
provides good performances in disturbances rejection and an efficient robustness to 
the process internal changes. The control of inverter-based microgrids has been deeply 
studied and some advanced methods have been successfully developed and tested (e.g. 
[3] [4] [5]). A preliminary work in power electronics [2] presents the successful applica- 
tion of the model- free control method to the control of dc/dc converters. This paper 
extends the previous results to the control of inverter-based microgrids in different 
situations relating to the islanded and grid-connected modes. In particular, we will 
show that the proposed control method is robust to strong load changes either in volt- 
age, current or power control. Simulations have been performed using the averaging 
method [6j with S$ - PE^\ [7J. 

The paper is structured as follows. Section II presents an overview of the model- 
free control methodology including its advantages in comparison with classical method- 
ologies. Section III discusses the application of the model-free control to inverters. 
Some concluding remarks may be found in Section IV. 

2 Model-free control: a brief overview 

2.1 General principles 
2.1.1 The ultra- local model 

We only assume that the plant behavior is well approximated in its operational range 
by a system of ordinary differential equations, which might be highly nonlinear and 
time- varying {*] The system, which is SISO, may be therefore described by the input- 
output equation 

E(t,y,y,...,yM,u,u,...,uM) = (1) 

where 

• u and y are the input and output variables, 

• E, which might be unknown, is assumed to be a sufficiently smooth function of 
its arguments. 

*Y i &—PES is a complete simulation platform that allows to perform co-simulation with C- written 
/ GNU Bison algorithms and SPICE OPUS; this software is particulary useful to simulate control 
law and power circuits whose topologies may switch during the simulation. 

tA new version that allows to drive SPICE OPUS under Matlab® has been released. 

*See [HE] for further details. 
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Assume that for some integer n, < n < l : ^ 0. From the implicit function 
theorem we may write locally 

= <£(*, y,y,..., j/"" 1 *, y (n+1) , . . . , y^,u, it,..., u^) 
By setting (£ = F + au we obtain the ultra-local model. 

Definition 2.1 /i/ //^ and y are respectively the variables of input and output of a 
system to be controlled, then this system admits the ultra-local model defined by: 

y( n) = F + au (2) 

where 

• a e K is a non-physical constant parameter, such that F and au are of the 
same magnitude; 

• the numerical value of F, which contains the whole ^structural information", is 
determined thanks to the knowledge ofu, a, and of the estimate of the derivative 

In all the numerous known examples it was possible to set n = 1 or 2. 
2.1.2 Numerical value of a 

Let us emphasize that one only needs to give an approximate numerical value to a. 
It would be meaningless to refer to a precise value of this parameter. 

2.2 Intelligent PI controllers 
2.2.1 Generalit ies 

Definition 2.2 [1] If n = 1, we close the loop via the intelligent PI controller, or 
i-PI controller, 

u = -- + V -+C(e) (3) 
a a 

where 

• is the output reference trajectory, which is determined via the rules of flatness- 
based control (73, HE/); 

• e = y* — y is the tracking error; 

• C(e) is of the form K P e + Kj J e. K P , Kj are the usual tuning gains. 

Equation is called model-free control law or model-free law. 

The i-PI controller [3] is compensating the poorly known term F. Controlling the 
system therefore boils down to the control of a precise and elementary pure integrator. 
The tuning of the gains K P and Kj becomes therefore quite straightforward. 
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2.2.2 Classic controllers 

See pj] for a comparison with classic PI controllers. 



2.2.3 Applicat ions 

See P2J [13l El [15l dU [171 [THl E] for already existing applications in various domains. 



2.3 Numerical differentiation of noisy signals 

Numerical differentiation, which is a classic field of investigation in engineering and 
in applied mathematics, is a key ingredient for implementing the feedback loop |3j 
Our solution has already played an important role in model-based nonlinear control 
and in signal processing (see [19] for further details and related references). 

The estimate of the 1 st order derivative of a noisy signal y reads (see, e.g., [20] ) 
y = — J Q T (T — 2t)y(t)dt where [0,T] is a quite "short" time window, § This window 
is sliding in order to get this estimate at each time instant. 

Denoising of y leads to the estimate y = ^ J Q (2T — 3t)y(t)dt 
The above results are the basis of our estimation techniques. Important theoretical 
developments, which are of utmost importance for the computer implementation, may 
be found in [21]. We establish the following hypothesis: 

Hypothesis 2.1 Jj$j The estimate of the derivative is realized with an high bandwidth 
and is not biased; in particular compared to the noise. 



2.4 A first academic example: a stable monovariable linear 
system 

Introduce as in [H [8] the stable transfer function 

(s + 2) 2 



(s + iy 



(4) 



2.4.1 A classic PID controller 



We apply the well known method due to Broi'da (see, e.g., [22] ) by approximating 
System [| via the following delay system 

Ke~ TS 



(Ts + 1) 



§It implies in other words that we obtain real-time techniques. 
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K = 4, T = 2.018, r = 0.2424 are obtained thanks to graphical techniques. The 
gain of the PID controller are then deduced [22]: K P = 1QQ 1 ^ Q ^ r ^ ) = 1.8181, Kj = 
1 - 0.7754, K D = e^r = 0.1766. 



1.33Kt w< ' ' J -^ u K 



2.4.2 i-PI. 

We are employing y = F + u and the i-PI controller 

u = -[F} e + f + C(s) 

where 

• i F ]e = [He ~ U, 

• y* is a reference trajectory, 

• £ = y* -y, 

• C(e) is an usual PI controller. 



2.4.3 Numerical simulations 

Figure 1(a) shows that the i-PI controller behaves only slightly better than the classic 
PID controller (Fig. 1(b)). When taking into account on the other hand the age- 
ing process and some fault accommodation there is a dramatic change of situation: 
Figure 1(c) indicates a clear cut superiority of our i-PI controller if the ageing pro- 
cess corresponds to a shift of the pole from 1 to 1.5, and if the previous graphical 
identification is not repeated (Fig. 1(d)). 



2.4.4 Some consequences 

• It might be useless to introduce delay systems of the type 

T(s)e~ Ls , TeR(s), L>0 (5) 

for tuning classic PID controllers, as often done today in spite of the quite 
involved identification procedure. 

• This example demonstrates also that the usual mathematical criteria for robust 
control become to a large irrelevant. 

• As also shown by this example some fault accommodation may also be achieved 
without having recourse to a general theory of diagnosis. 
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' 2 4 6 8 10 12 14 16 ' 5 10 15 

Time (s) Time (s) 

(a) i-PI control (b) PID control 




Time (s) Time (s) 



(c) i-PI control (d) PID control 

Figure 1: Stable linear monovariable system (Output (-); reference (- -); denoised 
output (. .)). 

3 Control of inverter-based microgrid 

We apply in this section the model-free control to the control of inverters used in typi- 
cal configurations within microgrid [23] in both stand-alone mode and grid-connected 
mode. We will show in simulation that the model-free control is robust to grid dis- 
turbances, load changes / load addition and efficient in direct abc frame (three-phase 
systems) control as well as in power control. 

All the inductors and capacitors described on the schemes have their values re- 
spectively close to 1 mH and 10 /xF. The dc bus voltage E is equal to 400 V. 

3.1 Voltage-controlled inverter 
3.1.1 Single load 

Consider a single-phase inverter working in stand-alone mode for which the output 
voltage is controlled (Fig. [2]). The load is a resistor R that switches from R « 10 Vt to 
R ?z 1000 at t = 0.02 s. Figure [3] presents the output voltage response V out of the 
inverter alone (without control - V* ut corresponds here to the PWM signal) and shows 
clearly that it exists a resonant effect after the load change (the LC filter becomes 
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less damped). Figures [i] and [5] present the output voltage response of the inverter 
according to the output voltage reference V* ut when a classical PI controller and an 
i-PI controller are considered. 




Figure 2: Full bridge inverter with a load. 




0.005 



0.015 
time (s) 



0.025 



0.03 



Figure 3: Response of the inverter alone. 
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Figure 4: Comparison between the PI and i-PI control for the voltage-controlled 
inverter (K p = 20, K { = 0). 
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It is difficult to adjust the PI controller in order to stabilize the system considering 
the load change. However, tuning the i-PI controller allows to adjust the dynamic 
transient response (Fig. [5]) especially during the starting transient of the inverter. 




Figure 5: Improvement of the transient dynamic performances using K p = 400, K{ = 
100. 



3.1.2 Multiple loads 

Consider a single-phase inverter working in stand-alone mode for which the output 
voltage is controlled (Fig. [6]). 




Figure 6: Full bridge with two loads. 



The inverter is firstly loaded by a resistor (load "1") and then a second unknown 
load (load "2") is added at t — 0.0042 s. Figure [7] shows the output voltage response 
of the inverter. Figure [8] presents the response for a PI control. Figure [9] presents the 
inverter output voltage response with an i-PI controller under different parameters. 
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Figure 7: Response of the inverter alone (the second load "2" is connected at t 
0.0042 s). 
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Figure 8: PI control (K p = 2). 
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Figure 9: i-PI control. 
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3.2 Three-phase current-controlled inverter 



Consider a three-phase inverter working in both stand-alone mode / grid-connected 



mode and controlled in current. The current in each phase is controlled (Fig. 10) by i 



*Z/2 a) -LL3CL 



PI (each phase has its own i-PI controller); I Lia 
going through the inductors L ia , L 2a , L 3a and 
reference currents. The load is composed of a three-phase resistor ( 
three-phase capacitor ( 



are the controlled currents 



^a^La^La are the corresponding 



10 Q) and a 



10 /jF). Figure 11 presents the voltages and currents of the 
inverter in stand-alone mode with a three-phase load change (R = 1000 C = 0.1/zF) 
at t — 0.012 s. Results are similar in the case of unbalanced conditions. 




3 Ph. grid 
3 Ph. switch 



Figure 10: Three-phase bridge inverter connected to the grid. 



A grid disconnection is presented Fig. 12 : a sinusoidal perturbation of 25 % of 



the grid amplitude at 500 Hz is added to the grid and the inverter is disconnected 
from the grid at t = 0.015 s. 
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Figure 11: Stand-alone mode controlled inverter with load changes. 
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Figure 12: Perturbated grid and disconnection from the grid. 
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3.3 Power-controlled inverter 



Controlling the output power of an inverter is important when considering paralleliza- 
tion of inverters and load sharing [21] [25] . 

Consider the single-phase full bridge inverter described Fig. [2] working in stand- 
alone mode; the load is a resistor (R = 100 Q). We consider in this section the control 
of the active power at the output of the inverter. The active power V is defined by : 



V 



J out 



(6) 



and its estimator is based on a moving- average filter. This is a direct control and 
the i-PI controller corrects the amplitude of the output sinusoidal signal in order to 
satisfy the power reference V*. Figure p~3| shows the active output estimated power of 
the inverter controlled by i-PI. A load change occurs (R = 50 Q) at t = 0.01 s. This 
strategy can also work in three-phase systems. 
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Figure 13: Power- controlled inverter. 



3.4 Parallel inverters 

Consider two single-phase inverters connected in parallel and working in stand-alone 



mode (Fig. 14). According to the power sharing methodology [26] . the inverter "1" 



is controlling the output voltage v out and the inverter "2" is controlling the current 
^L2- Figure [15] shows the output controlled voltage of the associated inverters. 
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Figure 14: Parallelization of inverters. 
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Figure 15: Output voltage of the parallelized inverters. 



14 



4 Concluding remarks 



We presented the model-free control methodology as a distributed control in an elec- 
trical network environment. Simulations show encouraging results and show that the 
model-free control has the following features : 

• robust to strong load changes (i.e. in case of the inverter loaded by a resistor, 
strong change of the resistor value); 

• robust to topological load changes (e.g. change of the load or addition of a load 
that may increase the order of the whole system); 

• robust to external perturbations (e.g. grid sinusoidal perturbation); 

• direct control in abc frame for three-phase systems; 

• non-linear control (e.g. power control). 

A combination of the proposed control strategies allows to extended the results to 
the control of multiple sources considering simultaneously voltage, current and power 
control. Further work will concern the study of the stability of the model-free control 
in networked systems. 
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